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his paper discusses the 
rotating stall in centrifugal 


T 


tional diagnostics and rotor dynamics 


compressors from vibra- 


standpoints. A diagnostic strategy is 
suggested from the overview of the 
case history and aerodynamic research. 
The literature survey shows insufficient 
understanding of rotor dynamic impli- 
cations of the rotating stall. Most of the 
literature on the subject addresses only 
aerodynamic aspects of the problem, 
and only a few papers treat a compressor 
as a coupled mechanical-aerodynamic 
system. The latter approach is exercised 
in this paper. The results of the 
rotor/fluid system Dynamic Stiffness 
identification are presented for the 
experimental centrifugal compressor, in 
rotating stall and in normal operating 
conditions. It is shown that the aero- 
dynamically-induced radial (direct 
spring) stiffness coefficients are 
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negative in normal operating conditions. 
The negative spring effect in the case 
of rotating stall is even more profound. 
The conclusion is that during rotating 
stall the rotor dynamic stability margin 
is reduced. 


Introduction 


The interest in radial compressor 
fluid-induced instabilities started in the 
mid-seventies. Since then, numerous 
case histories of high pressure, cen- 
trifugal compressor instabilities have 
been published [1-17]. Most of the 
publications report two types of rotor 
vibrational behavior: 1) High eccen- 
tricity, rotor first natural frequency re- 
excitation; and 2) Subsynchronous 
forward precession with rotative-speed- 
dependent frequency. The former is 
usually referred to as whip-type behavior 
and is normally associated with balance 
pistons, fluid-film bearings, and 


Impeller 









i XY 
Coupling Proximity 
probes 


Bronze 
bearing 


Figure 1. Experimental compressor. 
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labyrinth seals. The latter is called 
whirl-type behavior [27-28] and can 
be associated either with fluid-film 
bearings/seals or with rotating stall. 
Rotating stall usually results in an 
appearance of a low subsynchronous 
frequency component in the rotor 
vibration spectrum (frequency ratios 
are typically between 8 and 40% but 
can be as high as 80% of rotative 
speed, [20]). The emphasis on rotor 
behavior does not allow for details of 
particular flow patterns but treats the 
fluid (gas) “in average,” which falls in 
the scale of rotor motion. 


Experimental compressor 

An experimental rig (Figure 1 and 
Table 1) was constructed to evaluate 
the rotor/fluid system Dynamic 
Stiffness during rotating stall and 
normal operating conditions. The rig 
consisted of a single stage, centrifugal 
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Compressor type 

Working fluid 

Number of impeller blades 

Inlet diameter 

Outside diameter of impeller 
Diffuser/Impeller radius ratio 
Impeller tip/Diffuser radial clearance 
Impeller/Diffuser axial clearance 
Diffuser/Shroud clearance 


Discharge pipe dimensions 





Centrifugal, 1 stage, vaneless diffuser fully open, which allowed free flow 


Air at 84.1 kPa (12.2 psi), 20° C (68° F) 
ii 

51 mm (2 in) 

108 mm (4.3 in) 

1.8 

380 um (0.015 in) 

640 um (0.025 in) 

3 mm (0.12 in) 


with minimum restriction. A pressure 
transducer measured the compressor 
discharge pressure. 

A constant magnitude, rotating-force 
perturbator, used in modal testing, was 
attached to the nonrotating stub shaft. 
It transmitted an exciting rotating force 
with variable frequency to the compressor 
rotor. The perturbator consisted of an 


51 mm (2 in) inside diameter, 
360 mm (14.2 in) length 


outer rotating assembly supported by a 
rolling element bearing attached to the 


Table 1. Compressor design and operating parameters. 


compressor with a vaneless diffuser, 
driven by a 373 W (% hp) dc electric 
motor. The motor speed was controlled 
by a high output, variable power supply, 
and the speed could be controlled from 
near zero to 17,000 rpm. The motor 
was attached to the compressor rotor by 
a flexible disk coupling. 

The compressor rotor was separated 
into two pieces: a rotating section from 
the motor to the outboard end of the 
compressor, and a nonrotating stub 
shaft from the outboard compressor 
bearing to the outermost support bearing. 


The rotating section was supported 
by two bearings: a bronze bushing type 
at the inboard (motor) end, and a 
rolling element bearing at the outboard 
end of the compressor. The outer race 
of the bearing was installed in a collar, 
which was attached to the nonrotating 
outboard stub shaft. This allowed free 
rotation of the compressor rotor relative 
to the nonrotating stub shaft. A flexible 
membrane made of stiff rubber supported 
the stub shaft collar and the bearing 
assembly. Another bronze bearing 
supported the outboard end of the 
nonrotating stub shaft. 

An additional 0.8 kg (1.76 lb) disk 
was attached to the shaft between the 
inboard bearing and compressor to 
reduce the rotor system natural 
frequency. 


An orthogonal pair of eddy current 
displacement transducers was mounted 
between the air inlet and the impeller to 
measure the displacement of the rotor. 


A spring load frame was attached to 
the nonrotating stub shaft just outboard 
of the compressor. This was used to 
control the radial position of the 
compressor rotor and to counteract the 
gravity load of the rotor. 


A conical exhaust valve on a lead 
screw controlled airflow through the 
compressor. The valve was located at 
the end of the discharge pipe from the 
compressor (Figure 2). It could be 
adjusted through several turns from 
fully closed, which completely blocked 
the output flow of the compressor, to 


Diffuser 
scroll 


compressor support structure. Three 
radial springs with 120° angular spacing 
were attached from the outer rotating 
assembly to an inner rolling element 
bearing attached to the nonrotating stub 
shaft. One of the springs was tensioned 
more than the other two springs, and 
this created a net force in the direction 
of the most tensioned spring. As the 
perturbator outer assembly rotated, the 
springs rotated as well, and a rotating 
force was created. The force was 
transmitted through the inner rolling 
element bearing to the stub shaft, and 
then to the compressor rotor. A variable 
speed motor drove the perturbator via 
a flexible belt. The perturbator drive 
motor was equipped with a Keyphasor® 
probe to supply a phase reference and 
speed indication. The magnitude of the 
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Figure 2. The experimental compressor discharge flow/pressure control. 
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Figure 3. Experimental compressor performance chart. Q = rotative speed. 


rotating force, 1.3 N (0.3 Ib), did not 
change noticeably with speed, and the 
perturbator provided a variable speed, 
constant magnitude, rotating-force 
excitation to the compressor rotor. 

An experimental compressor 
performance chart is shown in Figure 3. 
It also depicts the points on the 
performance curves where the 
vibrational data was taken. 


The first series of experiments 
consisted of compressor startups with 
two different discharge valve positions, 
while the perturbator was at rest. Valve 
positions, which were as close to each 
other as possible, were chosen to 
provide normal operating conditions and 
rotating stall conditions. During normal 
operation, predominantly 1X forward 
vibration was observed, with no 


subsynchronous vibration. The data is 
presented in the form of a 1X, forward 
component Bode plot and full spectrum 
cascade plot (Figure 4). The Bode plot 
shows the existence of two rotor modes 
in the range from 0 to 17,000 rpm. 
With the help of an additional set of 
probes, the first mode was identified 
as an almost rigid body conical mode 
of the shaft, while the second was a 
bending mode. The data taken during 
rotating stall is shown in Figure 5. 
Since the Bode plot of the 1X forward 
component in this case does not differ 
strongly from the previous experiment, 
only a typical orbit and full spectrum 
cascade are presented. The experimen- 
tal compressor rotor exhibits rotating 
stall whirl with a subsynchronous fre- 
quency component of 0.21X in the for- 
ward direction. The subsynchronous 
component of the compressor rotor 
radial vibration is most likely caused 
by the diffuser rotating stall. 


Rotor dynamic characteristics of 
rotating stall 

The startup data does not provide 
sufficient information to identify any 
changes in the Dynamic Stiffness of 
the first mode due to rotating stall, 
because the mode occurs at rotative 
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Figure 4. Rotor lateral X and Y vibration data of the experimental compressor during startup without rotating stall: a) 1X forward 
component Bode plot, b) Full spectrum cascade plot. 
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Figure 5. Rotor lateral X and Y vibration data of the experimental compressor rotor during 
startup with rotating stall: a) Typical orbit/timebase, b) Full spectrum cascade. 


speed below the onset of rotating stall 
whirl. A series of nonsynchronous 
perturbation tests were performed at 
each of four compressor running 
speeds: 11,500; 13,500; 15,000; and 
17,000 rpm. At each running speed, 
perturbation trials were performed for 
the two valve positions used in the 
startup tests. The constant force pertur- 
bator generates a force of magnitude F 
with orientation ô rotating in the same 
direction as the compressor rotation 
(forward perturbation), resulting in the 
perturbation vector Fe”. The perturbator 
was ramped up from slow roll speed to 
3,000 rpm in order to identify Dynamic 


Stiffness of the rotor first mode. The 
vibration signals from two orthogonal 
eddy current displacement transducers, 
located at 90° to the right (X) and 0° (Y) 
relative to vertical (as seen from the 
driver), were filtered to the frequency œ 
of the perturbation force. This yielded 
amplitudes and phases for vertical and 
horizontal components of the response. 


In order to eliminate the minor 
influence of mechanical stiffness 
asymmetry, 1X vertical Y, and 
horizontal X o Vectors were combined 
into the circular component: 


Alfa) — Ý, = IX (1) 


oO 


Next, the perturbation process was 
repeated for perturbation in a direction 
opposite the rotation of the compressor 
rotor (reverse perturbation) for each 
running speed Q and valve setting. In 
this case, the forward (in the direction 
of perturbation) component of the 
filtered-to-perturbation-frequency 
response is opposite to the direction 
of rotation. 

The nonsynchronous Dynamic 
Stiffness was calculated by taking the 
ratio of the force input to the response 
vectors for each speed in the databases. 


Fe” (2) 

DS (œ) = qn 
DDS (@) = direct (DS) (3) 
QDS (œ) = quad (DS) (4) 


Since the total Dynamic Stiffness of 
the compressor consists of a mechanical 
and an aerodynamic part, the former 
was measured by running nonsynchro- 
nous perturbation on the compressor at 
zero rotative speed (DS,). The differ- 
ence between the Dynamic Stiffness 
of the running compressor and that of 
the stopped compressor represents 
aerodynamically-induced Dynamic 
Stiffness, DS ar 


DDS „(© )= direct( DS-DS) (5) 


air 


2 
= K air —M,;,;0 -2M Qo 


air2 


ODS „i, (©) = quad (DS — DS,) (6) 


= Dair -T 


air 


where Mgirp Kje and D ip 
inertia, radial stiffness, and damping, 
respectively, generated by the air/impeller 
interaction, Mg; is a fluidic mass term, 
and T represents a tangential force 


coefficient. 


are air 


Equations 5 and 6 show the parame- 
ters of Direct and Quadrature Stiffness 
extrapolation, that fit the data within an 
accuracy of about 2% compared to actual. 
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The raw data and the described data 
processing is illustrated in Figure 6. 
The results of all perturbation tests, 
reduced to the parameters of air-induced 
Dynamic Stiffness (Equations 5 and 6), 
are presented in Table 2. 


The most significant conclusion 
from the results in Table 2 is that aero- 
dynamic forces applied to the impeller 
generate negative spring stiffness 
(damping D,;, and stiffness K,;, 
appear negative in all experiments). It 
is also important that this stiffness 
drops even more during the rotating 
stall. This result is also shown in Figure 
7. Inconsistency in tangential force 
coefficient T suggests that it was too 
small in the experiment and was not 


determined with sufficient accuracy. 
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Case history 


A 5-stage propylene compressor 
provides an example of the industrial 
centrifugal compressor behavior under 
the conditions of surge/rotating stall. 
The orbit (Figure 8) shows a typical 
forward precession with frequency 
0.15X. This is supported by the lateral 
vibration, half spectrum cascade data 
(Figure 9). The gas pressure measured 
at the outlet of the third stage is pre- 
sented in the timebase waveform for 
the steady state regime (Figure 10) and 
half spectrum cascade plot for the 
startup (Figure 11). The steady state 
data shows the combination of vane 
passing high frequency with low 
frequency pulsations. Using the half 
spectrum cascade plot, the latter are 
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Figure 6. Identification of the air-induced Dynamic Stiffness component. 
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identified as a surge component (3.5 Hz) 
and rotating stall component (0.15X). 
It is important to note that the half 
spectrum cascade of the thrust probe 
data (Figure 12) shows a strong activity 
at the surge frequency, while the rotat- 
ing stall component is much smaller. 
The described behavior is a reflection 
of the surge, which propagates axially, 
while rotating stall is mainly active in 
the circumferential direction. 


Instrumentation and diagnostics 
strategy 

A typical machine train for a 
centrifugal compressor unit is shown in 
Figure 13. It includes an engine (in this 
case an aeroderivative gas turbine), 
gearbox, and the compressor itself. 
Different instrumentation configurations 
for vibration, as well as performance 
diagnostics and monitoring of the 
engine, are beyond the scope of this 
paper. For the centrifugal compressor, 
the instrumentation package normally 
includes a Keyphasor® transducer, two 
sets of proximity probes in XY config- 
uration installed at the inboard and out- 
board bearing locations (see 
Compressor DE and Compressor NDE 
in Figure 13), two axial transducers, 
and performance monitoring 
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Rotative speed Parameters 


Kair Ib/in 

Dair Ibes/in 

11,500 rpm Mair1 Ibes7/in 
Ibes’/in 

Ib/in 

i Ib/in 

i Ibes/in 

13,500 rpm Ibes’/in 
Ibes7/in 

Ib/in 


i Ib/in 

i Ibes/in 

15,000 rpm Ibes?/in 
Ibe°s?/in 

Ib/in 


Ib/in 
Ibes/in 
Ibs’/in 
Ibes7/in 


Ib/in 


17,000 rpm 


No 
rotating 
stall 


-52.16 
-3.4x10° 
-2.7x10° 
-9.75x10° 
2.83 
-134.6 
-7.0x10° 
-4.3x10° 
-10.0x10* 
-0.52 
-66.7 
-5.7x10° 
-4.3x10° 
-10.0x10° 
-0.52 
-63.2 
-5.2x10° 
-1.5x10° 
-8.7x10° 
3.6 


With 
rotating 
stall 


-63.2 
-3.0x10° 
-3.0x10°* 
-9.7x10° 
3.6 
-186.4 
-8.3x10° 
-7.7x10° 
-7.7x10° 
-0.07 
-134.1 
-5.7x10° 
-4.3x10° 
-10.0x10° 
-0.52 
-130.8 
-6.4x10° 
-5.0x10° 
-8.3x10° 
2.7 


Stiffness change, 
%, relative to 
normal conditions 





Table 2. Parameters of the aerodynamically-induced Dynamic Stiffness. 


instrumentation. As a minimum, the 
compressor performance monitoring 
instrumentation is supposed to include 
inlet and discharge pressures and 
temperatures, and discharge flow 
measurements. Using the described 
instrumentation, the rotating stall or 
surge can be identified based on the 
following information: 


e Radial vibration 
Steady state data: direct orbits, 
full spectrum 
Indication: subsynchronous 
forward component 


Possible causes: rotating stall, 


fluid-induced instability, and surge 


Transient data: direct orbits, 
full spectrum cascade 


Indication: subsynchronous 


forward component with frequency 


proportional to the rotative speed 
(0.1 to 0.8X) 


Possible causes: rotating stall, 


fluid-induced instability 


e Axial vibration 


Steady state data: spectrum 
Indication: subsynchronous 


component 


Possible causes: rotating stall 


(possible weak coupling), and 


surge 


Transient data: spectrum cascade 


Indication: subsynchronous 


component with low constant 


frequency 


Possible cause: surge 


e Pressure, temperature, flow 
Trend data (dc signal) 


Indication: pressure, head, and 
flow decrease, temperature 
increase 


Possible causes: surge, less likely 
rotating stall 


Note: if pressure head and flow 
indicate a position on the perform- 
ance curve close to the surge 
point, it could support rotating 
stall or surge versus fluid-induced 
instability 


Dynamic data 

Indication: vane passing frequency 
amplitude modulation of pressure, 
with constant flow 

Possible cause: rotating stall 


Indication: pressure vane passing 
frequency amplitude modulation 
with correlated flow oscillations 

Possible cause: surge 


As indicated above, a strong con- 
clusion can only be reached by the 
correlation of vibration data with 
performance information. It is important 
to note that rotating stall can originate 
in a particular stage of a multistage 
compressor. To localize the rotating 
stall source, additional pressure 
transducers may be required. 


Conclusions 


Although, the rotating stall aero- 
dynamic effects have been investigated 
extensively [16-23], the rotor dynamic 
implications of this phenomenon have 
received much less attention [24, 27- 
29]. A direct quotation from the paper 
[26] by Colding-Jorgensen serves as an 
example of insufficiency of rotor 
dynamic data: 

“In spite of all reservations due to 
difficulties and lack of data, it is the 
author’s opinion that the impeller- 
casing effect is mainly responsible for 
a negative direct stiffness effect.” 
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Figure 7. Aerodynamically-induced radial stiffness of experimental compressor. 
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Figure 9. Half spectrum cascade plot of the vertical inboard lateral vibrations of the rotor 
at the outboard bearing of the 5-stage propylene compressor during startup. 
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Figure 10. Gas pressure oscillations measured at the outlet of the third stage of 5-stage 
propylene compressor in steady state regime. Note high vane passing frequency 
modulated by the low frequency of rotating stall/surge. 
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Figure 8. Typical orbit of a 5-stage 
propylene centrifugal compressor rotor 
in rotating stall conditions. 


This paper provides some practical 


recommendations on the rotating stall 


diagnostics as well as some insight into 


rotor dynamic implications of the 


phenomenon. In short, the conclusions 


can be formulated as follows: 


The rotating stall conditions 
manifest in the rotor vibration 
signature as rotor forward pre- 
cession with subsynchronous fre- 
quency, which tracks the rotor 
speed. This behavior is referred 
to as a rotating stall whirl. 

From the diagnostic standpoint, 
rotating stall differs from the 
other instabilities of the whirl 
type in its strong dependence on 
the compressor operating condi- 
tions; it should disappear with 
increased flow. It differs from 
surge by its frequency propor- 
tionality to the rotative speed as 
well as by the predominantly 
axial character of surge-caused 
vibrations. 

The parameter clearly identified 
by nonsynchronous perturbation 
tests is aerodynamically-induced 
direct spring stiffness. In all 
conditions it was negative. From 
the practical standpoint, it means 
that even in normal operating 
conditions, centrifugal compressors 
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A native of St. Petersburg, Russia, Paul joined BRDRC in 
1990. He earned an M.S.M.E. and Ph.D. in Applied 
Mathematics from State Technical University of St. Petersburg. 
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